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Abstract 

We study a 2-spin quantum Turing architecture, in which discrete local rotations 
■ {cKm} of the Turing head spin alternate with quantum controlled NOT-operations. 

0> , We show that a single chaotic parameter input {am} leads to a chaotic dynamics 

I in the entire Hilbert space. The instability of periodic orbits on the Turing head 

^ ■ and 'chaos swapping' onto the Turing tape are demonstrated explicitly as well as 

, exponential parameter sensitivity of the Bures metric. 

o : 

^ : 1 Introduction 

^ ■ In recent years problems of quantum computing (QC) and information processing have 
Q\ . received increasing attention. To solve certain classes of problems in a potentially very 
O I powerful way, one tries to utilize in QC the quantum-mechanical superposition principle 
^ I and the (non-classical) entanglement an undertaking which, at the same time, should 
0\ • contribute to our basic understanding of quantum mechanics itself (see e.g. [0). However, 
building a large-scale quantum computer remains an extremely difficult task. The major 
obstacle is the coupling of the quantum computer to the environment, which tends to 
destroy quantum-mechanical superpositions very rapidly. This effect is usually referred to 
Ch I as decoherence. Present-day technology does not yet support the realization of a practical 
^ ' quantum computer. On the other hand, there might be interesting small-scale physics in a 
O^' pure quantum regime based on a few pseudo-spins (qubits), which are realizable right now. 
>• . Chaotic behaviour as an exponential sensitivity to initial conditions has been well es- 

^ , tablished in classical non-linear systems. The deterministic chaos, which occurs in non- 
^ I dissipative systems, can typicahy be found starting from regular states as a function of 
- - ' some external control parameter. However, there seems to be no direct analogue to chaos 
in the quantum world: If two quantum states are initially almost identical (that is, their 
scalar product is very close to 1), they will remain so forever, since the Hamiltonian evolu- 
tion is a unitary mapping which preserves scalar products. According to this negative result, 
the semi classical 'quantum chaology' has been constrained to studying some quantum- 
mechanical 'fingerprints of chaos' (like spectral properties), and non-trivial transitions from 
the quantum - to the classical domain and vice versa, following Bohr's correspondence prin- 
ciple (see e.g. [|, |^). In addition, experimental progress in mesoscopic physics, e.g. the 
transport of electrons through so-called 'chaotic quantum dots' 0, has attracted a great 
deal of interest, the results of which give numerical evidence for weak chaos (indicated by 
level repulsion) [|^, ^. 

While most models of QC have been concerned with networks of quantum gates, which 
are reminiscent of classical integrated circuits, models based on quantum Turing ma- 
chines (QTM) p, [ly] have been described along different lines but have not given rise 
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to much potential for future applications up to now. In both cases the complexity of the 
computation is characterized by sequences of unitary transformations (or the correspond- 
ing Hamiltonians H acting during finite time interval steps). The study of quantum chaos 
based on quantum gate networks has so far been proposed e.g. by implementing quantum 
baker's map on a 3-qubit NMR quantum computer |11[], by realizing a quantum- mechanical 



delta-kicked harmonic oscillator or a harmonically driven oscillator in an ion trap |12, 13 



and by showing quantum-mechanical localization of an ion in a trap fl^, . In all these 
cases some sort of sensitivity has been predicted with respect to parameters specifying the 
dynamics (e.g. the respective Hamiltonian) . 

Here we explicitly describe an iterative map with a few qubits which, though based on 
standard gates, can be thought to be realized as a QTM architecture: Local transformations 
of the Turing head controlled by a Fibonacci-like sequence of rotation angles alternate with 
a quantum-controlled NOT-operation with a second spin on the Turing tape. This chaotic 



control can generate a chaotic quantum propagation in the 'classical' regime |TH], defined 
here as one with the Turing head being restricted to an entanglement-free state sequence 
('primitive') [|T^. It will then be shown that chaos on the Turing head can be found also for 
the quantum-mechanical superposition of those primitives, implying entanglement between 
head and tape as a genuine quantum feature ('non-classical regime'). Finally, due to this 
quantum correlation, we find a chaotic propagation even in the reduced subspace of the 
Turing tape ('chaos swapping'), and as a result also in the total network state \il)n)- This 
behaviour should be contrasted with that of a regular QTM with its long-time revivals, 
which are absent in the 'chaotic' QTM. 



Description of chaotically driven quantum Turing 
machine 



The quantum network |T8[ to be considered is composed of (= M + 1) pseudo-spins 
Ip^(ai). p = 0, 1; /i = S', 1, 2, ■ ■ ■ , M (Turing-head 5*, Turing-tape spins 1, 2, ■ ■ ■ , M) so that its 
network state \ip) lives in the 2''^"'"^-dimensional Hilbert space spanned by the product wave- 
functions \j^^'^k^'^'> ■ ■ - l^^'^^) = \jk---l). Correspondingly, any (unitary) network operator 
can be expanded as a sum of product operators. The latter may be based on the following 
SU (2)-generators 

\ M _ pM I p(m) \ ip) _ ■ fbip) ■ nip) 

Ai — -TqI "t -'10 ? ^2 — ^-'01 "10 

\M _ pM _ p(m) I p(m) _ i(m) (t\ 

-^3 — Ai - -^00 ' ^0 - -^11 + -Hdo - -L 5 Uj 

where Pj^^'' = \pf'^^''^\q\ is a (local) transition operator. For simplicity we restrict ourselves 
here to M = 1. 

The initial state \ipo) will be taken to be a product of the Turing-head and tape wave- 
functions. For the discretized dynamical description of the QTM we identify the unitary 
operators t/„, n = 1,2,3, •• • (step number) with the local unitary transformation on the 
Turing-head S, Ui^, and the quantum-controlled-NOT (QCNOT) on {S, 1), U^^'^\ respec- 
tively, as follows: 

U2m^i = Ui^2 = i(^) COS {am/2) - xf^ % sin («^/2) (2) 
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f/(^^.i) = Pl^P Ai') + PIP = , (3) 



where the Turing head is externally driven by the Fibonacci-like sequence Om+i = «m+«m-i 
(mod 27r), ao = 0. The mth Fibonacci number is then controlled by ai via 

= - 7-) , (4) 

where P := , 7 := ■ The sequence {ttm} (mod 2tc) acts as a chaotic input 
(Lyapunov exponent: ln/5 > 0). It is useful for later calculations to note that 

= F(m) ■ /? + F(m - 1) , 7"^ = F(m) ■ 7 + F(m - 1) , (5) 

where F(m) := (/^"^ - 7™) /^S, the mth Fib onacci number with -F(l) = 1. The chaotic 
sequence of Fibonacci-type can be interpreted as a temporal random (chaotic) analogue to 
1-dimensional chaotic potentials in real space [0, [20[] . 

First, we restrict ourselves to the reduced state-space dynamics of the head S and 
tape-spin 1, respectively, 

Af V) = (^n|AP ® , A«(n) = (^„|i(^) A«|^n) . (6) 

\ipn) is the total network state at step n, X\^\n) are the respective Bloch- vectors. We intend 
to show that these local propagations are chaotic, too. Due to the entanglement between 
the head and tape, both will, in general, appear to be in a 'mixed-state', which means that 
the length of the Bloch- vectors in (|^) is less than 1. However, for specific initial states 
lipo) the state of head and tape will remain pure: As \^)^'^^ '■= (lO/^-* ± ll)''"'^'*) are the 

eigenstates of with A^^^jip) = ±\±f \ the QCNOT-operation U^^'^^ of equation (D 
cannot create any entanglement, irrespective of the head state \'^)^^\ i.e. 

f/(5-i)|^p)^ = xf\vf^®\-fK (7) 

As a consequence, the state \'il)n) remains a product state at any step n for the initial 
product states = \(pof^^ O |±)(^) with |(^oP = cos ((^o/2)|Op) - z sin ((^o/2)|lp) and 
the Turing head then performs a pure-state trajectory ('primitive') on the Bloch-circle 
{X['\n)=0) 

li^i) = I^^P ® \±f' , {Xi'\n\±)y + {Xi^\n\±)f = 1 . (8) 



Here A^ (^|±) denotes the Bloch- vector of the Turing head 5* conditioned by the initial 
state \il>o). From the Fibonacci relation and the property (0) it is found for Iv?^/"^^ ® 
\+f^\n = 2m, and ip^ = ao = that 

Af ^(2m|+) = sinC2^(+) , (2m|+) = -cosC2^(+) , (9) 
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where C2m{+) '■= ^^'^ ^^k\'^'^ ~ 1|+) — order to derive the cor- 

i=i 

responding expression of \^^\n\—) for If/?^/"^^ ® l"/^'*; "we utihze the following recursion 
relations for the cumulative rotation angle C„(— ) up to step n 

C2m{ — ) = —C2m-l{ — ) , C2m-l( — ) = + C2m~2{~) ■ (10) 

Then C2m{—),C2m-ii—), respectively, satisfy the following expressions: 

m 

C2U-) = -C2^-2(-)-a™ = (-1)™-'E(-I)'«i 

m 

C2m-lH = -C2^-3(-)+«,n = (-1)™E(-1)'«.' (H) 

yielding A2 (^|— ) = sinC„(— ), A3 {n\—) = — cosC„(— ) [cf. (p[)]. The Fibonacci property 
implies that both primitives, \(Pnf^'' ® l+Z^"* ^-^id Iv^^/'^^ ® l""^''; chaotically driven. 

From any initial state, lipo) = a^~^^\(fQ)^^'^ \+)^^^ + '^'■^■'Iv^o )'''^^ ® we then obtain 

at step n 

= a(+)|v9+/^) ® + a(-V;/^) ® 1-/1) (12) 

and, observing the orthogonality of the li/^-*, 

Af(n) = |aWpAf(n|+) + |a(-)rAf(n|-). (13) 

This trajectory of the Turing-head 5* represents a non-orthogonal pure-state decomposition. 
By using (D, (|Tl]), (0) (with a(+) = a^~'> = l/^/2) we finally have for |^o) = |0/^) ® |0p) 



(A2"^)(2m), A3'^''(2m)^ = cos^m ■ (sinSm, — cosS^) 
(Af)(2m - 1), Af)(2m - 1)) = cosS^ ■ (sin^^, -cos^^) , (14) 



where Am '■= oim + Oim-2 + • • • , '■= ttm-i + «m-3 + • • • • The cxprcssiou (|T4D indicates 



that also in the 'non-classical' regime the local dynamics of the Turing head is controlled 
by a 'chaotic' driving force, since the sequences Am and Bm, namely {a2m} (mod 27r) or 
{«2m-i} (mod 27r), respectively, are in fact both chaotic, as {am} (mod 27r) is. 

The Bloch- vector \^^\n) can alternatively be calculated directly from the initial state 
(here: |V^o) = |0p) ® \0p) and for any control angle Oi by using equation (|T^) and the 



relations 



A 



«L (^^m _ _ m = odd 



7 ' — yy m = even 



I ^ _ = even , ^ ' 

demonstrating a striking computational reducibility in that one needs for calculating A^'^-' (n) 
neither the total network state l^^^) nor to follow up each individual step n. 
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3 Instability of periodic orbits 



Now we verify that the periodic orbits on the plane |o, A2'^'', A3'^''| are unstable, which proves 
that the dynamics of the Turing head ('output') is indeed chaotic (figure |l|). Because of the 
alternating character of the dynamics, equations (0), (Rf), it suffices to check the periodicity 
only for step n = 2m: The periodic orbits for A2 (0) = 0, A3 (0) = —1 must obey two 

constraints, C2m{+) = C2m(— ) = 27rp, p G Z and Om+i = (mod 27r) for n = 2m + 1 (one 
concludes that ai must be a rational multiple of n). By using the Fibonacci numbers (|), 
one finds Cf^^ {+) in (P) and Cf^^ {—) in (pTf) , respectively, for period = 2m as 

Cf^'H = ^ + 1""-' + i-irV5) . (16) 

Then let us consider a small perturbation 6 of the initial phase angle ao = 0, implying 
l^oj-^^ = cos(5/2)|0) — zsin(5/2)|l) and a perturbed Fibonacci-like sequence {a'^} (mod 
2n): 

= 6, a'l = ai, a'2 = ai + 6, ■ ■ ■ . (17) 

Similarly to (|I6|), we obtain for this case C2„(±) = (±)-|- AC2m(±), where the deviation 
terms from the periodic orbits read, respectively. 

By using (|l^), we are able to represent the evolution of the perturbation, AX^^\n), at the 
n = 2m-th step for the Turing-head dynamics as 

( AXi'\2m) \ _ ( Mn \ ( AXi'\o) \ 
[ AAf (2m) )~[ M22)[ AA^ (0) ) ' 

where AAf^(O) = sin5, AAf^(O) = -cos5; AAf^(2m) = cos(5a„) sin ((5a^_i), AAf^(2m) 
= — cos(5a;m) cos((5am-i); ^11 = cos{6am) sm{6am-i) / sin 5, M22 = cos{6am) cos{6am-i) / 
cos 5, respectively. One easily shows [cf. equation (|^)] that 

lim Mil = F{m - 1) = ^ (/^'^^^ - 7™-^) lim M22 = 1 , (20) 

which means that Mn grows exponentially (note that (3 > 1, 17I < 1), and the periodic orbit 
on the Turing head is thus unstable to any small perturbation 6 in the external control. 
As an explicit example consider the case of ai = ^ with \iIjq) = \0f^^ ® |0/^): By using 

equation (|^) in equation (|1^) for C2^^{+) = 2ttp the periodic orbits with period = 2m have 
to obey F{m -|- 2) = 1 (mod 5). Together with the second condition, F{m + 1) = 1 (mod 
5) for n = 2m + 1, we obtain F(m — 1) = 1, F{m) = (mod 5). For C2^^(— ) = 27rp it 
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follows likewise that F{m — 1) = 1 (mod 5) for m = even, F{m — 1) = 4 (mod 5) for m = 
odd. Considering the common condition for both cases, F{m — 1) = 1 and F{m) = (mod 
5), we find the smallest m satisfying both conditions, namely m = 20, i.e. n = 2m = 40. 

liniMn = F(19) = 4181 > 1, see figure . 

Remarkably enough, the local dynamics of the Turing tape also contains some expo- 



nential sensitivity to initial conditions here: \ipo) = \0f^^ ® X^i\n 







cos — «! + ^Jr^^ n = 0, 1 (mod 4) 

n = 2, 3 (mod 4) 



(21) 



cos 



where ■= -^{P"^ ~ 7™); [a] := n, a = n + r, n G Z, < r < 1. One easily confirms 

that there is no periodic orbit with period 2m = 2 (mod 4). The following expression for 
the deviation from the periodic orbits with period 2m = (mod 4) at step n = 2m + 2 
holds: 



AA^^^(2m + 2) 



M-AA5^^(2), limM = ^ 



71 ^'^ 
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m+l 



sm a 



m+2) 



sm(ai 



(22) 



(2) = cos(q;i + 6) - cos(q;i), AA3^^(2m + 2) = cos(am+2 + 5m+i) - cos(am+2) 



with AA^^^ 

the perturbation S in (p!7|) does not enter into the Turing-tape dynamics until n = 2! 
Equation ( p2D shows the exponential instability of the periodic orbit. Indeed, the Turing 
tape exhibits chaos only by means of the entanglement with the head ('chaos swapping'), 
not result of a direct chaotic driving force. 



4 Exponential parameter sensitivity 

The distance between density operators, p and p', defined by the so-called Bures metric 

m 

D%> ■■= {{p - p'Y} ■ (23) 

lies, independent of the dimension of the Liouville space, between and 2 [the maximum 
(squared) distance of 2 applies to pure orthogonal states, = 2 (1 — |('?/'|'?/^')p) ]. This 
metric can be applied to the total-network-state space or any subspace. In any case it 
is a convenient additional means to characterize various QTMs: For = c^i (Lyapunov 
exponent = 0) [|l^ and any S the distance remains almost constant (see figure j^c?); for the 
Fibonacci-like driving force, on the other hand, (amip) = ctm, «m(p') = OLm{p) + (^m-^) 
obtain an initial exponential sensitivity, which is eventually constrained, though, by < 2 
(see figure ^ - c). It is instructive to consider another regular QTM controlled by the rule 
Om+i = 2Q;m — «m-i (Lyapuuov exponent = 0): am(p) = mai, am{p') = Oim{p) — (m — 1) 5. 
Here we observe a revival in the evolution of D^, which confirms that periodic orbits are 
stable (figure ^ci), whereas there is no revival for the above chaotic system. Finally we 
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display the evolution of for the total network state \ipn)i which also shows exponential 
sensitivity. The respective distances for tape-spin 1 are similar to those shown. 

The ultimate source of the present chaotic behaviour is that any small perturbation 5 
to the initial state, \ipo{5)), is connected with a perturbed unitary evolution, f/(5), which 
implies that the scalar product between different initial states (as a measure of distance) 
is no longer conserved under these evolutions: 

O' := I (^o(5) 1^(5) (J {mm)? , D' = 2{l-0'). (24) 

Thus the initial state is directly correlated to its unitary evolution, which can lead to the 
exponential sensitivity to initial condition, whereas there is no chaos in a generic quantum 
system evolving by a fixed tl even if characterized by chaotic input parameters. This O' 
reminds us immediately of the test function O = \{ip\V'' (t) U(t)\ip)\'^ [|^, where U,V are 
specified by slightly different external parameters: The corresponding parameter-sensitivity 
has been proposed as a measure to distinguish quantum chaos from regular quantum dy- 
namics. The origin of chaos in our QTM may thus be alternatively ascribed to a perturbed 
V = U{S) in the control (c/. the comment by R. Schack p3[|). 



5 Summary 

In conclusion, we have studied the quantum dynamics of a small chaotically driven QTM 
based on a decoherence-free Hamiltonian. Quantum chaos has been shown to occur as 
an exponential parameter-sensitivity and a cumulative loss of control in a pure quantum 
regime. This might be contrasted with the usual quantum chaology, which is concerned 
essentially with semiclassical spectrum analysis of classically chaotic systems (e.g. level 
spacing, spectral rigidity). As quantum features we utilized the superposition principle 
and the physics of entanglement. Our dynamical chaos manifests itself in the superposition 
and entanglement of a pair of 'classical' (i.e. unentangled) chaotic state-sequences. Due 
to the entanglement, we can see the chaos in any local Bloch-plane. This indicates that 
patterns in reduced Bloch-spheres (a quantum version of a Poincare-cut, figure |^) should 
be useful to characterize quantum chaos in a broad class of quantum networks: Here^ a 
periodic orbit would be represented by finite set of fixed points on the plane |o, A2'^'', A3'^ |. 
It is noteworthy that this kind of control loss is completely different from the typical 
control limit of a quantum network resulting from the exponential blow-up of Hilbert- 
space dimension in which the state evolves [|^. It is expected that a QTM architecture 



with a larger number of pseudo-spins on the Turing tape would also exhibit chaos under 
the same type of driving. However, it is just the chaos in small networks which might 
be interesting for experimental studies, especially in the form of ensembles thereof. The 
Fibonacci-like sequence should, however, be considered but a special example for chaotic 
input. Such inputs would, of course, have to be avoided in quantum computation; otherwise 
the resulting quantum dynamics would easily become chaotic in its entirety! 
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Figure |I|: A input-output scheme of our quantum Turing machine (QTM). 

Figure ^: Turing-head patterns {0, X2{n), ^sin)} under Fibonacci control for initial state 
l^o) = |0p) ® \0p. (a): ai = fyr (periodic); (6): ai = | x 3.141592654 (aperiodic) and 
total step number n = 10000. 

Figure |^: Evolution of the (squared) distance D'^^, between QTM state with (p') and 
without (p) perturbation 6. ai = §7?, {ipo) = \0f^'^ ® |0)(^) for p, and (^cos(5/2)|0)(^) - 



isin(5/2)|l/'^''j (g) \0f^^ for p'. (a): chaotic input according to equation (^ (inset shows 
initial behaviour in more detail; cf. Table |l]) for the Turing-head state, 6 = 0.001; (6): the 
same as (a) but for 6 = 0.00001; (c): the same as (a) for total network state \ipn)', (d): 
Dpp, within the Turing head subspace for am = ai (-D^ ~ 0, sohd line A, 6 = 0.001, 0.0005) 
and Q!m+i = ~ Q^m-i (dotted line B, 6 = 0.001; boxed line C, 6 = 0.0005), inset shows 
line B on larger scale. 
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Table 1: Analytical evolution of (squared) distance D^^, between the Turing-head states 
as of Fig l^a: With increasing time step n the originally small phase angle grows at a 
Fibonacci-like rate (separately for n = even and n = odd), such that shows exponential 
growth initially and then wild oscillations between and 2. 



n 







1 - cos(15) 


1 


1 - cos(15) 


2 


I + |cos(2ai + 26) - icos(2ai + 16) - icos(15) + |cos(2ai) 


3 


i - iCos(25) 


4 


i - iCos(25) 


5 


i - icos(2ai + 36) - icos(35) + icos(2ai + 26) 
— ■|cos(2ai — 16) — |cos(15) + ^cos(2a;i) 


6 


i + icos(6ai + 46) - |cos(6ai + 36) - icos(35) 
— ;|cos(6ai + 16) — icos(l(5) + |cos(6q;i) 


7 


1 — |cos(6ai + 56) — |cos(55) + |cos(6ai + 46) 
— icos(6ai — 16) — icos(15) + ^cos(6ai) 


8 


i + icos(8ai + 66) - icos(8ai + 55) - icos(55) 
— ;|cos(8ai + 15) — |cos(15) + ^cos(8a;i) 


9 


i - icos(8ai + 85) - |cos(85) + icos(8ai + 65) 
-icos(8ai - 25) - |cos(25) + icos(8ai) 


10 


i + icos(16ai + 105) - icos(16ai + 85) - icos(85) 
-icos(16ai + 25) - jCos(25) + icos(16ai) 


11 


i - icos(16ai + 135) - icos(135) + icos(16ai + 105) 
— |cos(16ai — 35) — ;|cos(35) + ^cos(16ai) 


12 


i + |cos(24ai + 165) - icos(24ai + 135) - icos(135) 
-icos(24ai + 35) - icos(35) + icos(24ai) 
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